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ABSTRACT 


Mean  flow  and  disturbance  measurements  have  been  made  on  a  flat 
plate  model,  one  meter  long  in  the  CWRU  low  speed  wind  tunnel  using 
hot-wire  anemometry.  Tests  were  conducted  on  a  smooth  plate  and  with 
two  different  grades  of  sandpaper  covering;  first  with  a  sandpaper 
covering  giving  Re^  values  centered  around  10  and  then  with  a  coarser 
sandpaper  having  Re^  values  of  about  150. 

The  mean  flow  measurements  indicated  Blasius  profiles  for  the 
smooth  plate  and  with  the  finer  roughness.  With  the  coarser  rough¬ 
ness,  the  profile  was  best  approximated  by  a  Blasius  profile  at  the 
early  stations  but  then  progressed  toward  a  turbulent  profile. 

The  finer  surface  roughness  (Rc^  1  10)  had  only  a  very  slight 

effect  on  the  u'  distributions  relative  to  those  for  a  smooth 
rms 

plate.  This  showed  up  primarily  through  examination  of  the  distri¬ 
bution  of  local,  turbulence  intensity,  u'  /u.  which  proved  to  be  a 

rms 

much  more  sensitive  indicator  than  dimensionless  turbulent  energy, 

2  2 

u'  /!'  .  for  the  coarser  roughness  (Re,  ~  150),  local  u'  intensi¬ 
ons  c  k  rms 

ties  in  the  boundary  layer  of  the  order  of  5%  were  observed  before 
the  mean  profile  departed  from  a  Blasius  shape.  The  progression 
toward  turbulent  mean  velocity  profiles  was  accompanied  by  large 
increases  in  turbulent  intensities  over  the  whole  boundary  layer  but 
especially  near  the  wall.  A  study  of  disturbance  spectra  for  this 


i  i 


<';isc  shows  largest  amplitudes  and  amplifications  at  frequencies  well 
below  the  Blasius  neutral  curve.  The  maximum  amplification  rates 
are  about  three  times  those  for  Tollmien-Schlicht infi  instability 
over  a  smooth  plate.  Transition  occurs  at  Re  '  350,000  which  is 

A 

much  below  that  for  a  smooth  surface. 
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CHAPTER  I 


INTRODUCTION 


On  a  smooth  surface,  transition  of  the  boundary  layer  presumably 
occurs  as  the  result  of  the  amplification  of  the  disturbances  ini¬ 
tially  present  in  the  flow.  In  this  model  of  the  transition  process, 
the  disturbances  amplified  are  those  whose  frequencies  lie  within 
the  neutral  curve  for  Tollmien-Schlichting  waves  and  are  dependent 
on  the  Reynolds  number  of  the  boundary  layer.  If  transition  on  a 

smooth  plate  occurs  at  distance  X  downstream  of  the  leading  edge 

co 

(with  corresponding  Reynolds  number.  Re  =  UX  /v),  then  for 

Co 

roughness,  transition  occurs  at  some  location  less  than  , 

o 

corresponding  to  a  transition  Reynolds  number,  Re^  =  UX^/v, which  is 

smaller  than  Re  (Fig.  1).  Roughness  enhances  friction  and  drag, 
o 

and  also  increases  the  momentum  thickness  of  the  boundary  layer. 

A  minimum  critical  roughness  is  sometimes  defined  as  one  for  which 

Re  is  approximately  0.95  Re  [1]. 

o 

The  treatment  of  roughness  is  divided  into  consideration  of 
single  roughness  elements,  both  two-dimensional  and  three-dimensional, 
and  of  distributed  roughness  which  is  considered  to  be  a  collection 
of  three-dimensional  protuberances  with  statistical  distribution  of 
height,  k,  and  distributed  over  the  plane  of  the  surface. 

There  are  many  experimental  results  for  single  roughness  ele¬ 
ments  and  there  is  even  a  basic  understanding  of  the  destabilizing 


1 


2 


mechanism.  On  the  other  hand,  much  less  has  been  determined  about 
distributed  roughness  and  there  is  no  understanding  of  its  destabi¬ 
lizing  mechanism. 

For  single  two-dimensional  roughness  elements,  Dryden  jl] 

presents  a  curve  of  Re  /Re  vs.  -rr  as  the  best  correlation  of  the 

t  t  <S* 

o 

body  of  data  to  that  time.  Smith  and  Clutter  [2]  argue  that  this 
is  too  simplistic  a  description  since  there  are  also  effects  of  the 
shape  of  the  element  and  the  length  of  laminar  flow  ahead  of  the 
element.  They  argue  that  a  more  fundamental  parameter  for  effects 

v 

in  the  vicinity  of  a  roughness  element  is  Re^  =  — .  They  found 

that  for  a  given  shape  of  roughness.  Re,  .  .  ,  was  found  to  be 

k, critical 

relatively  constant  and  independent  of  roughness  height  and  location, 
pressure  distribution  and  free-stream  turbulence  level.  There  was, 
however,  a  major  effect  of  roughness  shape. 

A  most  convincing  argument  on  the  mechanism  by  which  a  two- 
dimensional  roughness  element  destabilizes  a  laminar  boundary  layer 
is  by  Klebanoff  and  Tidstrom  [3].  They  show,  based  on  their  experi¬ 
ments,  that  the  basic  mechanism  by  which  a  roughness  element  induces 
earlier  transition  is  through  the  destabilizing  influence  of  the 
wake-like  recovery  zone  of  the  element  (Fig.  2).  Disturbances  in 
the  boundary  layer  are  greatly  amplified  in  the  recovery  zone.  The 
growth  rates  are  many  times  those  of  the  unaffected  Blasius  layer. 

The  fact  that  the  local  profiles  in  the  recovery  zone  depend  prima¬ 
rily  on  the  nature  of  the  roughness  element  establishes  Re^  as  the 
primary  parameter  of  the  phenomenon,  and  the  dependence  on  the  pro- 


perti.es  of  the  parent  flow  is  very  weak. 

For  three-dimensional  roughness  elements,  the  critical  values 
of  Re^  are  higher  than  for  two-dimensional  roughness  but  the  desta¬ 
bilizing  mechanism  may  be  similar  in  that  disturbances  can  be  great¬ 
ly  amplified  by  the  wake-like  unstable  profiles  downstream  of  the 
elements . 

In  view  of  the  results  for  single  roughness  elements  and  the 
known  sensitivity  of  boundary  layer  stability  to  small  profile 
changes,  the  analytic  investigations  to  date  of  distributed  roughness 
are  all  based  on  the  premise  that  distortion  of  the  mean  flow  by  the 
roughness  enhances  instability  and  leads  to  earlier  transition. 

Singh  and  Lunfley  (4  ]  made  calculations  to  find  the  roughness 
induced  mean  velocity  profile  which  they  expressed  in  terms  of  the 
spectral  density  of  the  roughness  and  an  influence  function.  The 
resulting  profile  had  an  inflection  point  deep  in  the  inner  viscous 
region  and  in  their  opinion  does  not  contribute  to  instability. 
Outside  the  wall  region,  the  profile  they  obtained  is  more  stable 
since  the  curvature  in  the  vicinity  of  the  critical  layer  has  been 
increased.  Since  experimentally,  roughness  leads  to  earlier  tran¬ 
sition,  Singh  and  Lumley  argue  that  even  if  the  profile  is  more 
stable,  the  roughness  also  introduces  disturbances  with  energy  at 
wavenumbers  to  which  the  profile  is  unstable. 

Lessen  and  Gangwani  [5]  model  the  problem  by  introducing  a  wavy 
wall  to  represent  one  component  of  the  roughness  spectrum.  The 
equation  for  this  roughness  induced  disturbance  is  the  Orr-Sommerfeld 


4 


equation  with  phase  velocity  equal  to  zero,  but  with  a  non-homoge- 
neous  boundary  condition  at  y  =  0.  They  retained  the  Reynolds 
stress  term  in  their  mean  flow  equation  and  so  calculated  the  dis¬ 
torted  mean  flow.  Once  this  distorted  profile  was  found,  the  stabi¬ 
lity  of  this  new  profile  was  examined  and  its  critical  Reynolds 
number  calculated.  The  new  profile  displayed  a  point  of  inflection 
near  the  wall.  The  new  profile  was  found  to  be  slightly  less  stable 
than  the  original  Blasius  profile.  The  minimum  critical  Reynolds 
number  is  reduced  by  about  52  in  the  example  presented  by  Lessen 
and  Gangwani.  This  amount  of  change  is  probably  not  very  signifi¬ 
cant  experimentally. 

Merkle,  Tzou  and  Kubota  [6]  assumed  that  roughness  enhances  mo¬ 
mentum  and  heat  transfer  near  the  surface  and  modeled  this  using  a 
turbulent  eddy  viscosity  in  the  vicinity  of  the  roughness,  in  a 
region  they  called  the  roughness  sublayer.  They  concluded  that  the 
roughness  sublayer  alter?  the  mean  velocity  and  temperature  profiles 
within  the  laminar  portion  of  the  boundary  layer.  They  obtained  an 
increase  in  shape  factor  primarily  through  increase  in  displacement 
thickness.  Stability  calculations  for  this  distorted  mean  profile 
indicated  reduced  stability.  An  attempt  by  the  authors  to  validate 
their  results  by  comparison  with  Achenbach's  experimental  results 
for  a  circular  cylinder  [7]  resulted  only  in  qualitative  agreement. 

Experimentally,  there  is  as  yet  no  demonstration  that  for  dis¬ 
tributed  roughness,  a  profile  change  is  responsible  for  enhanced 
instability.  The  present  work  intends  giving  a  definitive  experi- 


mental  look  at  this  question.  Accordingly,  the  mean  velocity  pro¬ 
files,  disturbance  energy  distributions  and  their  spectral  develop¬ 
ment  are  studied  in  order  to  help  understand  the  mechanism  by  which 
distributed  roughness  influences  transition.  The  investigation  is 
conducted  on  a  flat  plate  model  with  two  roughness  sizes  in  addition 


to  the  reference  case  of  a  smooth  surface. 


CHAPTER  II 


TEST  FACILITIES 

2.1  Low  Speed  Mind  Tunnel 

2.1.1  Physical  Descr i pt  ion 

The  low  speed  tunnel  is  of  the  open  circuit  type,  with  a  con¬ 
traction  ratio  of  10.4:1  (Fig.  3).  There  are  two  screens  ahead  of 
the  contraction,  each  with  an  openness  of  ft 5 % .  The  tunnel  is  ope¬ 
rated  by  a  53  IIP  DC  axial  fan.  The  test  section  is  7.5'  long  and 
its  cross-section  is  28"  x  28"  (Fig.  4).  Except  for  the  bottom  wall, 
the  rest  of  the  test  section  walls  are  made  of  transparent  plastic 
(plexiglas) . 

2.1.2  CjtJ  i b ra t  ion 

The  inside  of  the  tunnel  was  cleaned  and  sealed  to  achieve  as 
low  a  f roe-stream  turbulence  as  possible.  At  the  plate  location, 
values  of  free-stream  turbulence  intensity  below  0.03%  at  a  test 
speed  of  f>  m/s  or  were  obtained  as  shown  in  Figs.  5  to  7 .  The  hot¬ 
wire  readings  are  affected  by  the  velocity  components  which  are  per¬ 
pendicular  to  the  wire;  thus  rite  measured  turbulence  intensity  is: 


2 • 2  Elat  Plate  Model 
2.2.1  The  Plat e 

The  model  is  a  flat  smooth  plate,  41.5"  long,  25"  wide  and 
i/8"  thick  (Figs.  8-10).  The  leading  edge  has  a  1/32"  radius  fol- 
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lowed  by  a  smooth  taper  which  blends  into  the  flat  plate  surface 
1.2"  behind  the  leading  edge.  It  is  the  same  leading  edge  contour 
used  by  Strazisar  et  al.  [8,  9]  in  their  experimental  study  of  the 
stability  of  heated  laminar  boundary  layers  in  water.  The  surface 
of  the  plate  was  polished  to  a  mirror  finish  using  aluminum  polish 
compound  Met-All.  No.  1187. 

2.2.2  The  Mounting 

The  plate  is  supported  by  two  swept  front  supports  which  are  at 
a  70°  angle  to  the  plate,  and  two  rear  supports  perpendicular  to  the 
plate.  The  front  supports  have  a  symmetric  airfoil  cross-section 
with  2  3/4"  chord,  3/4"  thickness  and  the  cross-section  center  at 
the  plate  junction  is  16.5  cm  downstream  from  the  leading  edge  so 
as  to  prevent  separation  from  the  supports  and  generally  minimize 
creation  of  additional  disturbances  in  the  tunnel.  The  rear  sup¬ 
ports  are  3/4"  diameter  cylinders .  Their  center  is  80  cm  downstream 
of  the  leading  edge.  They  are  not  as  important  as  the  front  sup¬ 
ports  because  their  location  is  behind  the  test  region. 

The  choice  of  plate  location  was  based  on  the  wind  tunnel  ca¬ 
libration  results.  Since  the  longitudinal  variation  of  free-stream 
turbulence  intensity  is  constant,  the  leading  edge  of  the  plate  was 
located  near  the  door  of  the  test  section  for  convenience  in  making 
it  easier  to  clean  the  plate  surface.  It  was  positioned  at  a  17.5" 
height  above  the  floor  and  symmetrically  in  the  horizontal  plane. 

The  natural  vibrations  of  the  plate  on  its  supports  were  meas¬ 
ured  at  a  tunnel  speed  of  6  m/sec.  Fig.  11  shows  that  the  highest 
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amp  1 i tudo  appears  near  the  leading  edge  and  is  0.04  mm  and  occurs  at 
a  frequency  of  8  Hz  which  is  below  the  theoretical  critical  frequency 
hand  for  Tol Imien-Schlichting  waves.  Its  magnitude  is  much  below 
the  roughness  heights  tested.  The  amplitudes  at  the  middle  and  end 
of  the  plate  are  about  half  those  near  the  leading  edge.  A  flap  at 
a  1°  angle  into  the  test-stream  at  the  end  of  the  plate  was  used  to 
stabilize  the  stagnation  point  near  the  leading  edge,  in  order  to 
prevent  leading  edge  separation. 

2.2.1  Roughness 

Two  sizes  of  sandpaper  made  by  the  Norton  Company  were  used. 

They  were  applied  to  the  plate  by  using  RTV  adhesive.  The  finer 

sandpaper  (80  grit)  has  a  specified  nominal  average  grain  size  of 

266  microns  (0.0105  inch).  Thus  Re,  at  X  =  30  cm  at  U  =6  m /sec  is 

k  e 

10.6.  A  filler  was  applied  starting  1  cm  upstream  of  the  leading 
edge  of  the  sandpaper  and  blending  smoothly  into  the  sandpaper 
leading  edge. 

The  coarser  sandpaper  (24  grit)  had  a  nominal  average  grain 

size  of  1015  microns  (0.0408  inch)  so  that  Re,  at  X  =  30  cm  and 

k 

1)  =  6  m/sec  is  155.  In  this  case,  the  filler  began  1  inch  ahead  of 

the  leading  edge  of  the  sandpaper.  Because  of  the  different  thick¬ 
nesses  of  the  two  sandpapers,  the  ramp  angle  of  the  filler  was 
approximately  the  same  in  both  cases. 

Both  sandpapers  covered  the  whole  plate  surface  downstream  of 
their  leading  edges.  For  the  small  roughness,  Re^  was  calculated 


using  the  formula: 


Rek  = 


=  0.332  (— )  k  Xk 


(2.2) 


based  on  the  fact  that  for  small  heights  above  the  plate 
vX  -1/2 

— —  =  0.332  U  (— — )  iust  as  at  the  surface  [1].  Re  could  also 

dy  e  U  k 

R 

be  obtained  by  calculating  the  n  corresponding  to  k,  n  =  kw  — ,  and 

then  the  velocity  u,  could  be  taken  from  the  Blasius  curve.  This 
k 

method  was  applied  for  the  coarser  sandpaper  where  the  grains  are 
not  small. 

2 . 3  Measuring  Equipment 

2.3.1  Pitot-Static  Tube 

The  pitot-static  tube  used  was  of  the  Prandtl  type  made  by 
Sensor  United,  model  PAC-12KL  (Fig.  12).  It  is  inserted  from  the 
top  of  the  test  section  30  cm  downstream  from  the  leading  edge  and 
displaced  laterally  from  the  hot-wire  probe  by  5  cm  to  minimize  any 
interference  effect.  It  is  insensitive  to  yaw  or  pitch  angles  up 
to  30°  and  its  quoted  error  is  less  than  2%.  It  has  1/8"  outside 
diameter  and  is  positioned  5  cm  above  the  plate.  Thus,  the  distance 
from  the  surface  is  more  than  12  boundary  layer  thicknesses.  Since 
the  probe  is  also  more  than  10  probe  diameters  above  the  plate,  the 
effect  of  the  plate  on  the  readings  is  negligible. 

The  Reynolds  Number  based  on  the  impact  hole  diameter  (0.5  mm) 
for  6  m/sec  speed  is  approximately  100.  In  such  a  case  when  it  is 


larger  than  30,  viscous  effects  in  the  internal  flow  can  be  neglect¬ 
ed.  The  viscous  effects  in  the  external  flow  are  also  negligible 
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when  the  Reynolds  number  based  on  the  probe  radius  is  larger  than 
100,  which  is  the  case  here. 

For  T.I.  around  0.03%  and  dynamic  pressures  of  about  0.1  in. 

1  2 

H.jO.Lhe  Fluctuation  velocity  component  dynamic  pressure,  —  p(u') 

-9 

corresponds  approximately  to  9  x  10  in.l^O.  This  is  negligible 
relative  to  0.1  in.H^O.  Thus,  the  pitot  tube  reading  is  insensitive 
to  small  turbulence. 

The  time  constant  for  the  pitot  probe  is  approximately  five 
seconds . 

The  probe  is  connected  to  a  Meriam  manoaeter  which  it  ebl~ 
to  measure  0.00  1  in.il?0  (Fig.  15). 

2.1.2  Hot  W i re 

The  shape  of  the  hot-wire  probe  was  designed  appropriately  for 
boundary  layer  measurements  (Fig.  13).  It  consists  of  1/4"  stainless 
steel  tubing  in  which  electrical  wires  run  through,  that  are  con¬ 
nected  to  the  copper  prongs.  The  prongs  are  insulated  one  from  the 
other  by  5  minute  epoxy.  The  tips  of  the  prongs  are  connected  by 
Wollaston  wire  (90 7.  platinum,  10%  rhodium)  8  mm  long. 

The  length  of  the  control  portion  which  is  the  sensor  varies 
depending  on  the  etching  procedure  during  its  preparation  and  is 
approximately  0.8  mm  long.  Beneath  the  joining  of  the  prongs  there 
is  a  protruding  support  made  of  5  minute  epoxv.  From  cat hetometer 
measurements  supported  by  mean  velocity  profile  results,  the  wire 
is  0.2  mm  above  the  support  bottom. 

The  principle  of  the  hot-wire  and  the  factors  affecting  its 
operation  are  described  briefly  as  follows: 
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The  wire  is  heated  by  a  constant  current.  1 t  is  cooled  bv  the 
flow  causing  the  temperature  to  drop  and,consequently ,  the  electric 
resistance  to  diminish.  It  is  cooled  primarily  by  forced  convection 
and  heat  conduction.  Total  amount  of  heat  transferred  depends  on: 

1)  The  flow  velocity 

2)  The  difference  in  temperature  between  the  wire  and 
the  fluid 

3)  The  physical  properties  of  the  fluid 

4)  The  dimensions  and  physical  properties  of  the  wire 
The  last  three  are  known  and  thus  the  flow  velocity  can  be 

measured . 

The  hot-wire  anemometer  is  still  a  prime  instrument  for  tur¬ 
bulence  measurements.  Its  principal  advantage  is  the  separation  of 
the  mean  velocity  from  the  fluctuating  components  and  the  ability  Co 
measure  them  separately  at  the  same  time.  The  other  main  advantage 
is  in  its  being  thin.  The  Reynolds  number  based  on  wire  diameter 
is  1.34  for  a  speed  of  6  m/sec.  This  is  below  3,  which  is  the  limit 
for  smooth  flow  around  the  wire  without  separation  even  with  heating, 
because  Reynolds  number  decreases  with  temperature  |  JO]  (fig.  14).  The 
thermal  inertia  of  the  thin  wire  is  very  small  and  thus  the  response 
for  frequencies  of  interest  in  this  experiment  (up  to  5  KHz)  is 
practically  instantaneous.  It  catches  less  dirt  than  a  thick  sensor 
and  altogether  it  is  sufficiently  strong  and  rigid  to  sustain 
vibrations  and  stresses  caused  by  turbulence.  The  measured  velocity 
is  in  good  agreement  with  the  flow  velocity  as  measured  using  the 
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pitot-static  tube  for  turbulence  intensity  smaller  than  10%. 

The  joining  of  the  prongs  is  32,000  wire  diameters  behind  the 
wire  and  thus  there  is  no  real  influence  of  the  prongs  on  the  flow 
around  the  sensor  portion  of  the  wire. 

There  are  two  problems  witli  the  hot  wire.  The  first  is 
drifting  and  the  second  is  contamination.  Both  can  be  treated  prop¬ 
erly  as  will  be  explained  in  the  next  chapter.  More  about  hot  wires 
can  be  found  in  ref.  {11  1. 

The  hot-wire  anemometer  set  is  a  Shapiro  and  Edwards  model  50 
(Kig.  15).  It  includes  an  amplifier,  a  resistance  bridge,  a  poten¬ 
tiometer,  a  mean  square  output  meter  and  a  square  wave  generator. 

The  potentiometer  allows  measurement  of  wire  currents  up  to  10  ^  amp. 
and  so  enables  accurate  velocity  measurements.  A  compensation  cir¬ 
cuit  is  available  but  is  needed  only  for  frequencies  above  5KHz  and 
was  not  utilized  in  this  experiment. 

The  preparation  of  the  hot-wire  sensor  is  as  follows:  First, 
the  wire  is  soldered  to  the  prongs  and  then  the  probe  is  attached  to 
a  positioner.  The  wire  is  brought  into  contact  with  a  small  drop  of 
707.',  concentrated  nitric  acid  (pigs.  16-17).  An  electrolyzing  cur¬ 
rent  is  supplied  bv  a  1  1/2  volt  Eveready  battery  for  about  one 

minute.  Tiie  length  of  the  uncovered  (etched)  center  portion  of  the 
wire  is  a  function  of  the  drop  size,  and  is  of  the  order  of  0.8  mm. 
This  portion  which  is  the  sensor  has  a  diameter  of  0.00015  in. 

The  wire  is  then  checked  using  a  microscope  which  magnifies 
60  times.  It  is  important  to  tiave  a  straight  and  clean  wire. 
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The  soldering  of  the  wire  to  the  prongs  with  a  soft  solder  had 

to  be  done  quickly  in  order  not  to  create  excessive  heat  in  the  wire 

and  thus  avoid  any  elongation.  If  overheated,  the  uncovered  portion 
of  the  wire  could  not  sustain  as  much  stress  as  the  rest  of  the  wire, 

and  at  the  end  of  the  etching  process  would  become  deformed. 

The  other  factor  which  affects  wire  shape  is  the  surface  tension 
of  the  drop.  To  avoid  any  adverse  effect,  the  wire  has  to  be  taken 
out  of  the  drop  very  carefully. 

2.3.3  Traverse  Mechanism 

The  traverse  mechanism  (Fig.  13)  was  designed  to  give  accurate 
vertical  readings  of  the  probe  position.  It  has  a  metric  Mitutoyo 
micrometer  which  can  be  moved  longitudinally  with  the  help  of  a 
motor.  There  is  an  error  of  ^0.2  mm  in  longitudinal  positioning. 

The  micrometer  enables  vertical  precision  to  '0.01  mm  and  so  this  is 
the  accuracy  of  the  hot-wire  probe  position  because  it  is  attached 
to  tlie  micrometer.  The  height  can  be  read  on  the  micrometer  scale 
and  readjusted  to  zero  each  time  the  hot  wire  is  brought  to  the 
plate  surface.  The  precision  of  lateral  positioning  is  not  as  impor¬ 
tant  as  the  vertical  one  and  there  is  an  estimated  error  in  lateral 
position  of  at  most  0.5  mm. 

The  traverse  mechanism  is  outside  the  top  of  the  test  section 
and  the  hot-wire  probe  goes  through  and  moves  in  the  slot  made  in 
the  top  wall.  This  location  of  the  mechanism  was  chosen  because  it 
allows  the  probe  to  have  a  simple  shape,  has  a  shorter  shaft  than 
a  probe  from  the  side  wall  and  does  not  vibrate  at  the  test  speed. 
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(t  also  creates  less  disturbance  than  a  probe  from  the  side  wall. 

The  slot  was  placed  above  the  middle  of  the  plate  where  the 
turbulence  intensity  is  smallest. 

2.1.4  Spectrum  Analyzer 

A  Unigon  Model  4512  spectrum  analyzer  was  used  to  measure 
spectra.  It  has  an  averaging  capability  for  a  chosen  time  interval. 
The  bandwidth  and  the  amplitude  scale  could  be  adjusted.  A  Polaroid 
camera  attached  to  the  screen  enabled  pictures  of  the  spectrum 


to  be  taken. 


CHAPTER  III 


EXPERIMENTAL  PROCEDURE 


3.1  Ca  1  jLb  ration 

Calibration  of  the  hot  wire  was  performed  with  every  run 
because  the  cold  wire  resistance  was  changing  due  to  heat,  dust  and 
natural  deterioration.  To  eliminate  drifting,  it  was  necessary  to 
allow  two  hours  for  the  system  to  warm  up  before  using  the  hot  wire. 
The  contamination  problem  was  solved  by  changing  the  wire  every  few 
weeks . 

Calibration  was  always  performed  with  the  wire  at  3  cm  above 

the  plate,  and  30  cm  downstream  of  the  leading  edge.  The  cold  wire 

resistance  R  was  measured.  The  overheat  ratio  of  1.3  was  ad- 
o 

justed  by  setting  the  wire  resistance  R  to  a  value  of  1.3  R  . 

w  o 

The  wind  tunnel  speed  was  increased  gradually  and  each  time  the 

wire  current  and  the  manometer  readings  were  taken.  In  such  a 

manner,  the  calibration  curve  for  the  velocity  as  a  function  of  the 

wire  current  was  obtained  and  plotted  (Fig.  19).  It  was  almost  a 

straight  line  as  it  should  be  according  to  King's  law  [111: 

L2  =  a  +  b/u  (3.1) 

w 

After  the  calibration,  the  pitot  tube  was  pulled  out. 

3 . 2  Smooth  Plate 

3.2.1  Adjustment  of  Zero  Pressure  Gradient 

The  outer  flow  velocity  was  measured  along  the  plate  using  the 
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the  critical  stability  region;  that  is  where,  according  to  stability 
theory,  the  disturbances  are  amplified.  Since  tile  concern  was  to 
achieve  a  constant  pressure  in  this  region,  the  rear  supports  were 
shortened  by  I  cm,  and  so  the  plate  was  tilted  by  a  1°  angle.  The 
angle  of  attack  could  he  control  Led  further  by  elevating  the  rear  of 
the  plate  using  the  screws  in  the  rear  supports  but  this  was  not 
necossa  ry . 

As  can  be  seen  from  Fig.  22,  this  adjustment  yielded  a  constant 
pressure  region  between  X  =  30  cm  and  X  =  60  cm.  The  mean  velocity 
profiles  in  this  region  show  good  agreement  with  the  theoretical 
Kins  ins  curve  (Fig.  23).  In  this  region,  the  Reynolds  number  based 
on  displacement  thickness  is  390  to  83r'  which  is  in  the  critical 


stability  range  for  amplification  of  Tol lmien-Schl ichting  waves. 
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Thus,  the  region  between  X  =  30  cm  and  X  =  60  cm  was  chosen  as  the 
test  region. 

3.2.2  Measurement  of  Mean  and  Fluctuating  Information 

The  mean  and  the  fluctuating  information  was  taken  at  the  same 
time.  For  the  smooth  plate,  the  probe  was  lowered  until  the  support 
touched  the  plate.  This  was  repeated  several  times  and  the  micro¬ 
meter  readings  were  taken.  They  were  always  within  r0.03  mm  of  a 
mean  value  and  thus  could  be  considered  accurate  enough.  At  this 
position,  the  wire  was  approximately  0.2  mm  from  the  surface.  This 
was  the  closest  position  from  the  plate  and  there  the  first  readings 
were  taken.  Measurements  were  repeated  at  vertical  intervals  of 
0.3  mm  until  there  was  no  current  change,  indicating  the  outer  flow 
regime . 


At  each  point,  the  current,  the  series  resistance  and  the 
average  wire  voltage  V(e ' ) ^  were  measured.  This  enabled  calculation 
and  extraction  of  the  information  of  interest.  The  mean  velocity 
was  taken  from  the  calibration  curve  (Fig.  19)  for  the  corresponding 
current.  The  fluctuation  rms  velocity  component  was  calculated  as 
follows  [111: 


(u1) 


2  = 


U 


e ' )  ‘ 


(3.3) 


R  I  aEZ 
w  w 

where  S  = - is  a  hot-wire  sensitivity 

2u 


R  =  1.3  R  -  the  wire  resistance  at  overheat 
w  o 

I  - 
w 


the  wire  current 


a 


=  (K  -R  )/R  =  0.3 

w  o  o 

!■;  =  (1-.  )/(l+2'u)  where  r  =  R  /(R  +  R^)  and 

H  -  the  series  resistance 

1/(1+  /u  / fn)  where  H  is  the  corresponding  velocity 
o  o 

for  a  zero  current  and  extrapolated  from  the  cali¬ 
bration  curve 
,,  -  the  mean  velocitv 


The  volt.iye  '/( e ' )  “  was  measured  by  taking  approximately  50 


readings  at  5  second  intervals. 


By  calculating 


lurbance  energy 


;y  Vo*')2  £ 


n/=f 

V(u')  =  u  ,  it  was  easy  to  get  the  dis- 

rms 


and  the  free-stream  turbulence  intensity 


V<V)2/u 


or  the  local  intensity 


V(u’)2/u. 


1.  1  Rougli  Plate 

1.  1.1  Additional  Adjustments 

The  finer  sandpaper  was  at  first  positioned  to  begin  at  3  cm 
downstream  of  the  leading  edge.  The  Re^  and  k/A*  at  this  location 
are  33  and  0.56  respectively.  The  boundary  layer  was  apparently 
tripped  by  the  leading  edge  of  the  sandpaper  and  the  flow  at  the 
l  irst  measuring  station, X  =  30  cm , was  already  turbulent.  The  paper 
was  then  cut  back  to  6.5  cm  downstream  of  the  leading  edge  of  the 
plate.  The  corresponding  Re^  and  k/j’*  were  23  and  0.38  respecti¬ 
vely.  Here  the  flow  was  laminar  at  the  first  station  and  for  this 
location  of  the  leading  edge  of  the  sandpaper,  the  measurements 
were  taken. 

The  coarser  sandpaper  began  initially  at  6.5  cm  from  the  lead- 
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ing  edge  of  the  plate  but  again  the  boundary  layer  was  tripped. 

The  corresponding  Re^  and  k/6*  were  314  and  1.48  respectively. 

The  flow  was  turbulent  even  at  X  =  20  cm.  When  the  paper  was  cut 
back  to  18.3  cm,  the  flow  at  the  first  station  was  laminar.  At 
X  =  18.3  cm,  Re^  =  192  and  k/S*  =  0.88.  This  was  the  location  of 
the  leading  edge  of  the  coarser  sandpaper  for  the  measurements 
presented  herein. 

3 . 3.2  Measurements  Taken  with  Each  of  the  Two  Sandpapers 

The  measurements  and  the  procedure  for  the  finer  sandpaper  was 
similar  to  those  for  the  smooth  plate.  For  measurements  with 
coarser  sandpaper,  precautions  had  to  be  taken  to  prevent  wire 
damage . 

A  glass  plate  of  1.4  mm  thickness  was  placed  on  the  surface 
under  the  probe  and  the  probe  was  lowered  until  the  support  touched 
the  glass  plate.  The  micrometer  reading  was  taken.  The  probe  was 
raised  and  the  glass  plate  removed.  The  probe  was  lowered  again  by 
1.4  mm  beneath  the  reading.  At  all  stations,  it  was  possible  to 
lower  the  probe  an  additional  0.5  mm  without  damaging  the  wire.  The 
probe  location  was  illuminated,  and  it  was  always  possible  to  see  a 
gap  between  the  wire  and  the  grit.  This  was  checked  many  times.  This 
is  possible  because  of  the  non-uniform  distribution  of  the  grains 
in  the  sandpaper  matrix  and  also  because  of  some  grains  which  were 
positioned  above  the  average.  If  the  plate  rested  on  those  higher 
grains,  then  the  wire  could  be  lowered  beneath  their  height.  We 
could  assume  that  there  were  enough  grains  of  such  size  to  allow 
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tii  is  to  occur.  Tit  us ,  the  first  reading  for  the  roughest  sandpaper 
was  taken  with  the  wire  at  0.3  mm  below  the  highest  grains  in  the 
near  vicinity  where  the  measurements  were  taken. 

In  addition  to  the  measurements  as  taken  earlier,  special 
data  were  taken  for  the  coarser  sandpaper.  Here  the  information 
was  taken  also  at  X  =  20  cm.  The  hot-wire  probe  was  positioned  at 
the  peak  energy  location.  The  amplitudes  were  measured  from  Pola¬ 
roid  pictures  of  the  analyzer  output  and  the  amplification  rates 
were  calculated  from: 


a . 

l 


1  dA 
A  dX 


(3.4) 

The  amplitude  gradient  was  measured  from  the  amplitude  curves. 

Data  were  also  taken  using  the  filters  on  the  hot-wire  anemo¬ 
meter  panel.  At  the  peak  energy  location  for  each  station,  the 
energy  was  measured  for  different  frequency  bands.  The  filters  were 


1-10,  10-100,  100-5000  Hz. 


CHAPTER  IV 


RESULTS  AND  DISCUSSION 


4 . 1  Smooth  P la t e 

To  provide  a  reference  for  the  later  measurements  with  surface 
roughness,  a  set  of  mean  and  fluctuating  measurements  was  taken 
on  the  smooth  aluminum  plate.  The  mean  velocity  profiles  at  X  =  30, 
40,  10  and  60  cm  downstream  of  the  leading  edge  are  shown  in 
ligs.  24  (a-d) .  Two  data  sets  taken  at  different  times  are  shown. 

The  Blasius  curve  taken  from  Sehlichting  [12]  is  plotted  for  com¬ 
parison  on  each  of  the  figures.  The  points  agree  with  the  Blasius 
curve  within  the  experimental  error.  The  profiles  corroborate 
the  earlier  indication  from  surface  pressure  measurements  that  the 
region  30  •  x  60  cm  is  a  constant  pressure  region. 

One  set  of  test  results  for  the  fluctuation  energy  nondimension- 

2 

a  I i zed  with  respect  to  U~  versus  dimensional  boundary  layer  coordi¬ 
nate  v  and  nnndimons ional  n  is  presented  in  Tigs.  25a  and  25b 
res pec t i ve 1 v . 

The  free-stream  turbulence  and  the  wind  tunnel  fan  rotation 
affect  the  magnitude  of  the  disturbance  energy,  but  the  shape  and 
order  of  magnitude  of  the  curves  remains  consistent.  Even  though 
the  data  points  are  not  always  on  the  curve  lines,  the  deviation  is 
within  the  experimental  error  and  the  curve  lines  are  definitely 
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rharaeterist ic  of  the  phenomena  taking  place  according  to  other 
test  results  taken  in  the  present  investigation.  The  shape  is  si¬ 
milar  to  that  obtained  in  other  experimental  studies  for  a  smooth 
plate,  for  example  that  of  Strazisar  et  al.  [8,9]. 

The  integrals  under  the  curves  and  the  peak  values  are  increas¬ 
ing  consistently  and  gradually  in  the  downstream  direction.  The 
energy  content  is  distributed  over  the  whole  boundary  layer  but  is 
concentrated  mainly  in  the  middle.  It  is  almost  symmetrically  dis¬ 
tributed.  The  peak  moves  away  from  the  wall  but  remains  at 

n  2.3  or  y  /<*>  ~  0.4  3  and  y  /6  *1.37.  The  peak  energy  is 
max  max 

increasing  approximately  linearly  with  distance  downstream. 

The  local  disturbance  intensity  which  is  with  respect  to  the  lo- 
ca I  velocity  is  plotted  in  Tig.  26.  It  is  seen  to  decrease  near  the 
wall  witli  downstream  distance.  The  energy  integrals  and  peak  values 
are  increasing  but  very  little.  Both  the  energy  and  the  local 
intensity  are  increasing  in  the  vicinity  of  the  edge  of  the  boundary 
layer,  thus  indicating  that  the  response  to  free-stream  disturbances 
is  not  confined  to  tiie  interior  of  the  boundary  layer.  The  turbu¬ 
lence  intensity  at  the  outer  edge  is  seen  to  increase  about  0.02% 
with  every  10  cm  downstream  from  the  leading  edge. 

The  reasons  for  showing  both  energy  and  local  turbulence  inten¬ 
sity  will  become  apparent  later  in  the  cases  with  roughness. 

h .2  Distributed  Roughness  of  Re^  (X  =  30  cm)  =  10,6 

Similar  measurements  were  taken  for  this  roughness  as  for  the 
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smooth  surface.  The  mean  velocity  profiles  at  X  =  30,  40,  50  and 
60  cm  downstream  of  the  leading  edge  are  shown  in  Figs.  27  (a-d) . 
dust  as  for  the  smooth  case,  a  number  of  data  sets  were  taken  and 
the  Bias i us  curve  is  plotted  for  comparison  on  each  of  the  figures. 
The  zero  reference  height,  as  for  the  smooth  wall  case,  was  assumed 
to  be  0.2  mm  beneath  the  hot-wire  when  the  probe  support  touched 
the  surface.  With  this  choice  of  n  =  0  (y  =  0) ,  the  Blasius  curve 
fits  the  data  points  within  the  experimental  error. 

One  set  of  test  results  for  the  energy  is  presented  in 
Figs.  28  a,b.  The  shape  and  order  of  magnitude  of  the  curves  for  this 
run  remains  consistent  for  the  results  of  other  tests  taken  in  this 
investigation  but  not  shown.  The  energy  has  the  same  shape  and  the 
same  gradual  increase  as  for  the  smooth  case.  It  is  again  mostly 
concentrated  at  the  center  of  the  boundary  layer  and  the  peak 
remains  at  n  '  2.3  as  for  the  smooth  surface. 

Wo  can  see  already  slightly  more  curvature  near  the  wall  than 
in  the  smooth  wall  case  (Fig.  25b)  as  the  flow  proceeds  in  the 
downstream  direction.  The  energy  at  the  outer  edge  of  the  boundary 
layer  is  also  slightly  larger  than  for  the  smooth  plate.  Within  the 
boundary  layer  the  energy  levels  are  slightly  smaller  than  for  the 
smooth  surface  perhaps  because  the  speed  was  lower  by  0.2  m/sec  and 
possibly  because  the  flow  was  less  turbulent  initially,  thus  affecting 
the  initial  energy  levels.  Although  the  energy  peak  at  X  =30  cm  is 
smaller  than  for  the  smooth  wall  case,  the  energy  growth  -  the 
ratio  of  energy  at  the  last  station  to  that  at  the  first  station  - 


is  107.  larger  than  before  and  that ,  perhaps,  is  an  indication  of  a 
slight  roughness  effect. 

A  more  pronounced  indication  of  the  roughness  effect  is  shown 
by  the  local  turbulence  intensity  distribution  presented  in  Fig.  29. 
Even  though  the  energy  was  smaller,  the  intensities  were  larger  than 
those  for  the  smooth  wall  (Fig.  26)  in  the  vicinity  of  the  wall. 

The  intensities  of  the  peaks  close  to  the  wall  are  increasing  with 
downstream  distance  in  contrast  to  the  results  for  the  smooth 
surface . 

In  summary,  this  small  roughness  (beginning  from  Re^  =  22.8 
and  k/(S*  =  0.38  at  X  =  6.5  cm)  has  only  a  very  small  effect.  It 
doesn't  perturb  the  flow  to  the  nonlinear  regime  and  both  the  mean 
velocity  and  the  disturbance  energy  profiles  resemble  those  of 
laminar  flow  on  a  smooth  plate.  Although  the  overall  effects  are 
small,  we  could  detect  changes  in  the  local  turbulence  intensity 
distributions,  particularly  in  the  vicinity  of  the  wall. 

4.  1  Distributed  Roughness  of  Re^  (X  =  30  cm)  =  155 

4.3.1  Mean  and  Fluctuation  Measurements 

A  set  of  mean  velocity  profiles  is  shown  in  Figs.  30  (a  -  d) . 

In  this  case  the  roughness  is  large  enough  so  that  the  effective 
zero  velocity  location  for  the  boundary  layer  could  not  be  deter¬ 
mined  a  priori.  The  data  point  closest  to  the  wall  was  0.3  mm 
below  the  Lop  of  the  largest  grains.  This  point  was  chosen  as  the 
reference  location  for  the  rest  of  the  data  taken.  The  nominal 
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gri  t size  is  known  and  thus  the  distance  of  the  reference  location 
point  from  the  nominal  bottom  of  the  grains  is  also  known.  By 
fitting  a  curve  through  the  data  points,  the  extrapolated  zero  velo¬ 
city  point  can  be  located.  For  the  firs't  two  stations  (X  =  30  cm, 

X  =  40  cm),  a  Blasius  curve  could  reasonably  be  put  through  the 
data  points  and  the  extrapolated  zero  velocity  locations  were  at 
approximately  the  same  distance  from  the  nominal  bottom  of  the  grit 
(0.2  mm).  An  attempt  to  fit  the  data  points  with  a  similarity  pro¬ 
file  for  8  =  -0.1  was  unsuccessful.  These  two  measured  profiles  are 
much  closer  to  the  Blasius  curve  and  compare  very  well  with  the 
measured  data  on  the  smooth  plate  once  an  adjustment  is  made  for 
the  extrapolated  zero. 

Curves  through  the  mean  velocity  data  points  for  the  two  down¬ 
stream  stations  (X  =  50  cm,  X  =  00  cm)  are  shown  in  Figs.  30c  and 
30 d .  The  extrapolated  zero  velocity  location  is  chosen  to  be  the 
same  as  in  Figs.  30a  and  30b.  Blasius  curves  are  drawn  for  compa¬ 
rison.  At  X  =  50  cm  the  mean  velocity  profile  has  changed  signifi¬ 
cantly  from  that  at  X  =  40  cm.  At  X  =  60  cm  there  is  the  beginning 
of  the  emergence  of  a  turbulent  profile  shape. 

if  there  are  inflection  points  near  the  wall  in  Figs.  30a  and 
30b,  it  is  practically  impossible  to  detect  them  considering  the 
experimental  error.  Tt  is  also  very  difficult  to  Jo  a  detailed 
search  in  the  vicinity  of  the  grit  because  of  prospective  wire 
damage  and  probe  interference  effects. 
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Fluctuation  data  for  two  separate  tests  taken  a  week  apart  are 
shown  in  Figs.  31a,  b.  The  first  test  results  are  indicated  by  I 
and  the  second  by  II.  The  first  test  corresponds  to  the  mean  velo¬ 
city  results  shown  earlier  (Fig.  30).*  The  n  =  0  locations  used  in 
these  figures  for  each  station  were  taken  from  the  extrapolated  zero 
velocity  points  of  the  mean  velocity  profiles.  The  shape  of  the 
energy  distributions  is  not  smooth  but  somewhat  oscillatory.  Both 
tests  show  more  than  one  local  peak  and  the  number  increases  with 
the  distance  downstream.  The  first  two  stations  display  more  or  less 
their  previous  shape.  Their  peaks  are  in  the  vicinity  of  n  *  2.3 
and  only  the  second  test  deviates  from  this  slightlv  at  the  second 
station.  The  downstream  profiles  here  are  much  fuller  than  for  the 
smooth  plate  or  the  finer  roughness  and  there  is  significant 
increase  in  energy  content  near  the  wall. 

For  the  finer  roughness  the  only  clear  indication  of  a  rough¬ 
ness  effect  was  through  the  changes  in  the  turbulence  intensity  dis¬ 
tributions.  Here  we  see  a  large  effect  on  the  energy  distributions 
as  well.  The  peaks  are  moving  closer  to  the  wall  and  most  of  the 
energy  content  is  concentrated  in  the  n  <  4  region  for  the  first 
three  stations.  Farther  downstream  there  is  also  significant 
increase  in  energy  content  at  the  outer  edge  of  the  boundary  layer. 


*  Mean  velocity  profiles  for  test  II  (4/21/80)  were  measured  at 
X  -  30  cm  and  X  =  40  cm.  they  are  verv  close  to  the  profiles  of 
Figures  30a  and  30b .  The  profiles  at  X  =  50  cm  and  X  =  60  cm  differ 
slightlv  from  those  of  Figures  30c  and  30d  in  a  manner  consistent 
with  the  different  progress i ons  to  transition  on  those  two  days. 


27 


In  between  X  =  40  cm  and  X  ~  50  cm  whore  a  velocity  profile  change 
takes  place,  the  energy  increases  by  an  order  of  magnitude.  Figs. 
32a, b  present  the  local  turbulence  intensity  for  these  two  tests 
respectively,  again  using  the  extrapolated  zero  from  Test  1.  The 
local  intensity  rises  to  more  than  15%  between  X=A0cm  and  X  =  50  cm 
emphasizing  that  the  disturbances  are  in  the  non-linear  regime. 

In  Test  I,  there  is  not  much  overall  energy  change  between  the 
last  two  stations  but  there  is  a  noticeable  buildup  in  the  outer 
portion  of  the  boundary  layer.  Accompanying  this  is  a  definite 
change  in  the  mean  velocity  profile  (Figs.  30c, d)  and  an  increase 
in  the  boundary  layer  thickness  even  though  it  is  not  yet  a  fully 
turbulent  flow. 

Test  II  shows  a  flow  that  seems  less  disturbed  initially.  Less 
of  an  increase  in  energy  occurs  between  X  =  AO  cm  and  X  =  50  cm. 

As  a  consequence,  there  is  still  an  increase  in  energy  between  the 
the  last  two  stations.  Only  the  last  station  indicates  a  signifi¬ 
cant  energy  increase  in  the  outer  portion  of  the  boundary  layer. 

The  mean  velocity  profiles  which  are  not  presented  here,  again 
indicate  a  Blasius  profile  for  the  first  two  stations  and  a  very 
slight  change  relative  to  Test  f  at  the  third  station.  The  inten¬ 
sity  peak  at  the  last  station  reaches  the  same  magnitude  as  at  the 
last  two  stations  of  Test  I. 

The  intensity  shape  is  bumpy.  However,  the  main  peaks  are  approach¬ 
ing  the  wall  t  consistent  with  the  observation  for  the  Re^  =  10.6  case. 
The  peaks  are  within  the  region  n  •-  0.5.  At  the  measurement  point 
closest  to  the  wall,  the  total  amplitude  ratio  relative  to  free- 


stream  turbulence  level  (assumed  to  be  0.02%)  for  both  experiments 
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is  already  about  1500  which  is  c  .  The  boundary  layer  at  this 
point  is  however  not  yet  fully  turbulent.  Maybe  the  ratio  is  larger 
closer  to  the  wall  hut  it  was  impossible  to  measure  any  closer. 

The  consistency  of  the  peak  intensity  development  shown  by 
Figs.  32a,  b  with  that  of  the  finer  roughness  indicates  that  the 
extrapolated  zero  point  is  reasonably  chosen.  The  development  of 
the  energy  and  the  intensity  indicate  that  the  fully  turbulent 
regime  is  not  yet  reached. 

An  attempt  to  isolate  a  roughness  effect  for  Test  II  by 
subtracting  Fig.  26  (the  smooth  plate  intensity)  from  Fig.  32b  (the 
rough  plate  intensity)  is  presented  in  Fig.  33.  The  curves  are 
essentially  the  same  as  in  Fig.  32b  indicating  that  the  overwhelming 
contributions  to  the  data  in  Fig.  32b  are  due  to  the  roughness. 

Fig.  34  shows  the  variations  of  the  peak  local  intensity  and 
the  square  root  of  the  peak  nondimcnsional  energy  with  Re^.  These 
results  are  consistent  with  what  was  said  earlier:  (a)  the  intensity 
is  more  sensitive  than  the  energy  and  thus  a  better  indication  of 
the  degree  to  which  the  flow  is  disturbed;  (b)  the  flow  in  Test  I 
Is  more  disturbed. 

In  summary,  the  intensity  magnitudes  for  both  experiments  are 
above  15%  at  the  third  station  (X  =  50  cm)  and  thus  are  clearly  non¬ 
linear.  The  intensity  for  the  second  station  is  4.5  -  8%  and  thus 
significant  too  and  probably  non-linear.  The  main  conclusion  is 
that  the  disturbances  grow  to  non-linear  amplitudes  before  detect¬ 
able  velocity  profiles  are  observed. 


4.3.2  The  Spectrum 

In  addition  to  mean  and  fluctuation  measurements,  one  set 
of  tests  was  made  using  the  llnigon  spectrum  analyzer.  Figs.  33  -  39 
show  the  spectra  taken  at  X  =  20,  30,  40,  30  and  60  cm  with  the 
coarsest  paper  and  Li  ^  =  6  m/sec.  For  X  =  20  cm  to  X  =  40  cm 
(Figs.  35,  36,  37)  gain  settings  were  inadequate  to  get  data  above 
20  Hz  from  those  pictures.  It  was  possible  to  read  tile  amplitude 
above  20  Hz  for  the  last  two  stations  but  the  accuracy  is  doubt¬ 
ful,  and  since  the  amplitude  slope  is  very  sensitive  to  any  small 
amplitude  change,  the  calculated  v  values  are  doubtful  too. 
Nevertheless  it  is  possible  to  extract  some  important  information 
from  this  test. 

Amplitude  distributions  for  frequencies  up  to  17.2  Hz  are  shown 
in  Fig.  40.  They  have  a  shape  resembling  what  is  obtained  in  linear 
stability  studies.  Fig.  41  presents  the  calculated  amplification 
rates  For  17.2  Hz  amplification  occurs  at  Re  .  -  330.  As  the 

i  6  * 

Re  increases,  successively  lower  frequencies  become  amplified  but 
their  maximum  amplification  rate  is  diminishing.  The  highest  ampli¬ 
fication  rate  for  17.2  Hz  is  located  at  Re  ,  =  650.  Lower  frequen- 

<5* 

cies  show  peak  rt.  at  successively  larger  values  of  Re  .  These 
1  (S  * 

peaks  occur  in  the  non-linear  region  of  the  flow  as  identified  from 

energy  and  intensity  distributions.  Fig.  42  shows  again  that  as 

Re  ^  inc revises ,  the  lower  i requeue ies  become  important.  A  plot  of 

“j  vrsus  Re.A  (Fig.  43>  indicates  lhat  higher  frequencies  (above 
max 
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17.2  Hz)  may  be  important  below  Re  =  650.  Note  that  the  Tollmien- 

o* 

Schlichting  hand  at  Re  .  =  650  is  between  45  -  85  Hz.  In  fact  if 

o  * 

we  extrapolate  the  trends  of  Fig.  41  to  higher  frequencies,  it  would 
suggest  the  higher  frequencies  become  more  important  close  to  the 
leading  edge. 

Ampl  if  icat  ion  rates  -r*.  for  the  frequencies  up  to  130  Hz  at 
X  =  50  cm  are  shown  in  Fig.  44.  Similar  behavior  of  a.  is  found  at 

i 

X  =  60  cm.  Amplification  is  shown  at  all  frequencies.  There  does 
not  seem  to  be  anything  special  about  the  Tol lmien-Schlichting 
hand.  Again,  at  X  =  50  cm,  the  disturbances  are  non-linear. 

The  experimental  "neutral  curve"  for  frequencies  up  to  17.2  Hz 
is  below  the  neutral  curve  for  the  linear  stability  of  non-parallel 
flat  plate  flow  as  calculated  hy  Saric  f l 3 ]  (Fig.  45).  It  is  also 
below  the  neutral  curve  for  the  linear  stability  of  non-parallel 
I’a  1  kner-Skan  flow  corresponding  to  6  =  -0.1.  This  corroborates  the 
earlier  indication  that  disturbance  growth  was  governed  by  non¬ 
linear  growth  effects  rather  than  by  linear  instability  of  an  in¬ 
flectional  mean  velocity  profile  which  hasn't  been  detected  in  any 
case . 

From  the  trend  of  the  curve  we  can  definitely  say  that  there  is 

earlier  amplification  and  very  likely  earlier  transition  than  for 

the  smooth  surface.  While  the  linear  theory  for  a  smooth  plate 

predicts  :l.  for  Re  ^  above  684,  here  it  is  larger  much  below  this 
ma  x 

number.  The  order  of  magnitude  of  (a.  6*  )  is  also  totally 

1  r  r  ma  x 


d i f  ferent  . 
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-a.  6* 

The  maximum  value  of  (  —  yt)  for  non-parallel  smooth  flat 

66  -6 

plate  flow  is  approximately  9.55  x  10  and  occurs  at  Re^  =  795  [9]. 
For  the  same  Re the  present  data  for  the  rough  case  shows  it  equal 
2.67  x  10  ^  and  the  value  of  is  even  higher  for  lower  Re^.  This 
means  faster  transition.  Considering  the  state  of  the  boundary 


layer  at  X  =  60  cm  (Re^  =  835),  transition  might  be  completed  at 

Re  '  350,000  which  is  very  low  compared  to  the  value  of  3  x  10^ 

X 

that  might  be  expected  on  a  smooth  plate  for  the  low  free-stream 
turbulence  intensity  of  the  tunnel. 

Frequencies  corresponding  to  roughness  wavelengths  are  of  order 


—  =  (~~) w  •  ln  the  present  experiment  such  frequencies  would  be  in 
the  range  of  1.5  to  j  KHz.  No  growth  was  observed  in  these  frequency 
bands . 


4.3.3  Fi 1  to  ring 

Another  experiment  was  performed  using  the  built  in  filters  on 
the  hot-wire  anemometer  set  at  the  peak  energy  location  at  X  =  30, 
40,  50  and  60  cm  downstream  from  the  leading  edge.  The  filters  were 
l  —  1 0 ,  10-100,  100-5000  Hz.  Fig.  46  shows  the  energy  results  for 


these  frequency  levels. 

There  is  a  shift  in  energy  to  higher  frequencies.  Most  of  the 
energy  in  the  non-linear  regime  is  concentrated  in  the  100-5000  Hz 
band.  This  didn't  show  up  from  the  spectrum  because  even  though 
the  amplitude  in  the  high  frequency  range  is  small,  its  integral  is 
significant  and  very  sensitive  to  a  small  amplitude  increase,  since 
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the  band  width  of  the  100-3000  Hz  band  is  so  large.  Major  high 
frequency  content  is  probably  in  the  100-200  Hz  band.  The  mechanism 
is  driven  by  frequencies  up  to  about  140  Hz  as  indicated  by  the 
spectra  (Figs.  3d,  39),  but  the  energy  is  dissipated  into  higher 
f  requenc ies . 

4 . 4  Leading  Edge  of  the  Dist ributcd  Roughness 

As  described  in  the  preceding  chapter  for  both  roughness  sizes, 
the  leading  edge  of  the  sandpaper  had  to  be  cut  back  in  order  that 
the  boundary  layer  not  be  tripped  by  the  leading  edge  of  the  sand¬ 
paper.  Such  tripping  did  occur  before  cutting  back  the  sandpaper, 
despite  the  use  of  a  filler  to  avoid  a  step  at  the  leading  edge  of 
the  sandpaper. 

The  finer  roughness  sandpaper  was  cut  from  X  =  3  cm  with  the 
corresponding  Re  =  33  and  k/<S*  =  0.56,  to  X  =  6.5  cm  downstream  from 

K 

the  leading  edge  and  the  corresponding  Re^  =  22.8  and  k/6*  =  0.38. 
Before  the  cutting  hack  the  sandpaper,  the  flow  was  turbulent  over 
the  whole  test  region  and  after  cutting  it  back  it  was  laminar. 

The  coarser  roughness  paper  was  cut  back  from  X  =  6  cm  with  the 
corresponding  Ke^  =  314  and  k/<S*  =  1.48, to  X  =  18.3  cm  with  corres¬ 
ponding  Re  =  198  and  k/fi*  -  0.88.  Before  cutting  the  sandpaper  back, 
transition  was  detected  already  at  X  =  20  cm  and  turbulent  flow  oc¬ 
curred  in  the  whole  test  region.  After  cutting  back  the  sandpaper, 
laminar  mean  profiles  were  obtained  at  least  to  X  =  40  cm  (Fig.  30). 


Eiven  though  Re^  was  equal  to  198  at  the  leading  edge  of  the 
coarser  sandpaper  and  only  33  at  the  leading  edge  of  the  finer  sand¬ 
paper  and  the  filler  ramp  angle  was  approximately  the  same  in  both 
cases,  the  latter  tripped  the  boundary  layer  while  the  former  did 
not.  Thus  Ke  is  not  the  only  parameter  affecting  the  flow.  The 

K 

same  can  probably  be  said  about  k/o*.  Thus,  the  response  of  the 
flow  to  the  leading  edge  of  a  rougness  region  cannot  be  simply 
described  in  terms  of  these  two  parameters.  The  filler  strip 
geometry  and  the  smoothness  achieved  are  important  factors  in 
influencing  the  flow. 
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CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 

A  preliminary  experimental  study  has  been  completed  of  the  dis¬ 
turbances  in  a  laminar  boundary  layer  due  to  distributed  surface 
roughness.  Tests  were  conducted  on  a  smooth  plate  and  with  two  dif¬ 
ferent  grades  of  sandpaper  covering. 

Over  the  portion  of  the  smooth  plate  displaying  zero  pressure 

/ 

gradient  (30  cm  <  x  <  60  cm),  the  mean  flow  velocity  profiles  agree 

with  the  Blasius  profile  within  the  experimental  error.  The  u' 

rms 

distributions  display  a  peak  at  n  *  2.3  that  increases  slightly  with 

distance  downstream  reaching  a  level  of  0jOO7  Ug  at  X  =  60  cm. 

For  distributed  roughness  having  Re^  ®  10,  the  results  are  quite 
comparable  to  those  for  the  smooth  plate.  The  mean  flow  velocity 
profiles  agree  with  the  Blasius  curve  within  the  experimental  error 
and  the  u^.^  distributions  are  closely  comparable  in  level  and 
distribution  to  those  for  the  smooth  plate.  There  is  however  indi¬ 
cation  of  a  slight  increase  in  local  intensity  in  the  very  near  vici¬ 
nity  of  the  wall.  It  is  well  to  point  out  that  although  the  dimen- 

2 

sionless  turbulent  energy  (u'  /U  )  and  the  local  turbulent  inten- 

rms  e 

sity  (u,rmg/u)  involve  the  same  fluctuation  amplitudes,  the  local 
turbulent  intensity  Is  a  much  more  sensitive  indicator  of  effects  of 
roughness. 
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For  distributed  roughness  with  Re  is  150,  the  picture  is  dif- 

K. 

ferent.  At  X  =  30  cm  and  at  X  =  40  cm,  the  mean  velocity  profile  is 
best  fit  by  a  Blasius  profile  whose  effective  origin  is  somewhere 
within  the  roughness  height.  Beyond  X  =  40  cm  the  mean  profile 
seems  to  progress  toward  a  turbulent  profile.  Up  to  X  =  40  cm,  the 
u'  still  peaks  at  n  ;  2.3  but  the  local  intensity  is  now  of  the 
order  of  5%.  Beyond  X  =  40  cm  there  is  a  large  increase  in  the  u' 
levels  and  an  increase  in  the  number  of  peaks.  Very  large  local 
intensities  are  obtained  in  the  near  vicinity  of  the  wall  (n  <  0.5) 
but  there  is  also  noticeable  buildup  over  the  entire  boundary  layer 
to  the  outer  edge.  A  study  of  disturbance  spectra  for  this  case 
shows  largest  amplitudes  and  amplifications  at  frequencies  well 
below  the  Blasius  neutral  curve.  The  maximum  amplification  rates 
are  about  three  times  as  large  as  those  for  Tollmien-Schlichting 
instability  over  a  smooth  plate.  Transition  occurs  at 
RCj.  '•  350.000  which  is  much  below  the  value  of  3  x  10^  or  so  that 
might  be  expected  for  a  smooth  surface  at  the  turbulence  level  of 
the  tunnel. 

To  conclude,  one  observes  that  disturbances  can  grow  to  non¬ 
linear  magnitudes  before  velocity  profile  changes  are  observed. 

The  indicated  importance  of  frequencies  below  those  of  the  Blasius 
neutral  curve  was  not  expected  and  requires  more  careful  examination. 
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Sdicm.it  ic  diagram  of  effect  of  distributed  roughness  on 
transition  on  a  f]at  plate. 


The  mean  velocity  distortion  due  to  2-D  roughness 
el emont . 


Tlu*  front  viow  of  tlu*  model. 
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Figure  14.  Flow  Regimes  About  Unheated  Circular  Cy¬ 
linders  of  High  Leng th- to-D lame  ter  Ratio 
Versus  Reynolds  Number  (Morkovin,  1964). 
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f-'ig.  24c.  Mean  velocity  profile  for  a  smooth  plate,  X  =  50  cm. 


8^-2 


59 


r  i 

S: 

V)  , 


I 

1  <o 

■f  - 

'  '  VQ  k) 


Ha  i 

^  pv 


;  Ci ;  <o 1 
j  cvy 


+  i  O 


<0 


V 


(t) 


I 

I 

I 

I 


^  ^  fO ' 


II 

o 


o 

nj 


<r, 

o 


U-4 

C 

u 

CL 


>, 

U 

o 

c 


oc 

•r4 

iu 


Y 


63 


£  o 

A  -<\  o 

<sj'  'a- 04 
-r  i  .  m  • 

<s*  n  ^ 

K  •  <j 


* 

«£ 

ti 

VJ 

s 

«»»• 

-J 

0) 


b<s> 


Fig.  27b.  Mean  velocity  profile  for  a  rough  plate,  Re  (X  =  30  cm)  =  10.6;  X  =  40  cm. 
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30  cm) 


profile  for  a  rough  plate.  Re.  (X  =  30  cm)  =  153;  X  =  40  cm 


Tost  I  -  Loc.il  turbulence  intensity  profiles  for  a  rough  plate, 
Ke  (X  =  TO  cm)  =  135. 
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li'.sL  (I  -  boc.il  turbulence  intensity  due  to  roughness  only, 
Kok  (X  =  30  cm)  =  15r>. 
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The  liif;lu'st  amplification  rate  versus  X. 
Ue,  (X  =  30  cm)  =  155. 


Arapl i f lea t i on 
Re,  (X  =  30  cm 


Non-para  1  le  1  theoretic.il  neutral-stability 
curve  for  (  =-0.1,  Saric  (1975) 

Non-parallel  theoretical  neutral-stability 
diagram  for  (•  =  0,  Saric  (1975) 

Roughness  results  for  (■  =  0,  Re^  (X  =  30  c 


600  700  soo  900  #e 

Neutral  stability  results. 


